Precisely timing dissipative quantum information processing 
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Dissipative engineering constitutes a framework within which quantum information processing protocols are 
powered by system-environment interaction rather than by unitary dynamics alone. This framework embraces 
noise as a resource, and consequently, offers a number of advantages compared to one based on unitary dynamics 
alone, e.g., that the protocols are typically independent of the initial state of the system. However, the time in- 
dependent nature of this scheme makes it difficult to imagine precisely timed sequential operations, conditional 
measurements or error correction. In this work, we provide a path around these challenges, by introducing basic 
dissipative gadgets which allow us to precisely initiate, trigger and time dissipative operations, while keeping 
the system Liouvillian time-independent. These gadgets open up novel perspectives for thinking of timed dis- 
sipative quantum information processing. As an example, we sketch how measurement based computation can 
be simulated in the dissipative setting. 



One of the main goals in quantum information science and 
in quantum technologies in general is to understand the infor- 
mation processing power available within the framework of 
quantum mechanics. Several quantum computational models 
have been proposed and analyzed, and each has its strengths 
and weaknesses. Every such model will suffer severely from 
the problem of decoherence due to noise from a practically 
inevitable environment, so that a lot of effort has been in- 
vested in developing methods for isolating systems from their 
destructive environment. In such a picture, quantum noise 
and dissipation is seen as necessarily being detrimental for 
coherent quantum state manipulation. Recently, however, a 
new avenue has been suggested for addressing this issue: ac- 
tively exploiting the dissipation into the environment. In a 
number of recent studies, it has been shown, both experimen- 
tally |T] |2 and theoretically OlfTTll that several of the basic 
tasks in quantum information science (quantum computation 
0, state preparation 0HS1 [Zl [HI , distillation (9), storage ifTOlL 
and cooling |[TT1l ) can be performed by engineering the sys- 
tem environment coupling, instead of isolating a subsystem 
and performing coherent controlled unitary dynamics on it. 

The most counterintuitive application of such a paradigm 
might be dissipative quantum computing (DQC) 0; a model 
of quantum computation in which one assumes a system- 
environment interaction described by a Markovian master 
equation, where the computation is encoded in the Lindblad 
operators and the outcome of the computation is encoded in 
the unique stationary state of the open system. The main 
benefit with this approach (that under appropriate conditions 
the unique stationary state is reached rapidly starting from 
any initial state) appears to imply a potentially severe disad- 
vantage: that operations cannot be performed sequentially in 
time. Steps of a procedure cannot be conditioned on previous 
steps, it is not clear how to "stop preparing" a state, and it 
is far from clear how to incorporate error correction into any 
such scheme. 

In this work, we open up new perspectives for dissipative 
quantum information processing by introducing and analyz- 



ing a number of dissipative "gadgets". They can be combined 
to act as time triggers for various dissipative operations; i.e., 
they allow a certain dissipative process to start (or stop) acting 
after a very specific point in time: They serve as convenient 
"clocks" in a framework with time-independent Liouvillians. 
The functioning of these tools is based on a counterintuitive 
property of classical Markov chains called the "cutoff phe- 
nomenon" lfT2l . which can be lifted up into the quantum set- 
ting [14]. We also show how such gadgets can be used in 
order to "translate" any scheme involving unitary dynamics, 
measurements and conditional dynamics into a time indepen- 
dent dissipative setting. We then outline how these time trig- 
gers can be used to perform dissipative measurement based 
quantum computation lfl3l . 

Although universal quantum computation can be performed 
dissipatively without recourse to timing of operations, or ini- 
tializations, we show here that using these simple and intuitive 
gadgets, one can obtain results which may be difficult to ob- 
tain otherwise. Indeed, the timer gadgets can be used to show 
that DQC is universal also for geometrically local interactions 
(6-local if embedded in a 3D lattice). In the original proof 
of DQC 0, the interactions were made 5-local, but in a way 
which cannot be embedded in a lattice geometry. 

We show that timed dissipation-driven quantum informa- 
tion processing is indeed possible. Needless to say, there is 
a lot of room to reduce the significant overhead involved in 
order to make the schemes more suitable for experiments. 

The paper is set up as follows. We first introduce an ini- 
tialization gadget which allows to prepare a given state for 
a finite amount of time. Next, we combine the initialization 
gadget with a special behavior from Markov chain mixing to 
construct a rudimentary timer gadget, which allows to trigger 
a dissipative operation at a specific point in time. Rigorous er- 
ror estimates for these two gadgets are given in the Appendix. 
We go on to outline an application of these gadgets: a dissi- 
pative one-way quantum computation (D1WQC) scheme. 

Basic building blocks. The basic conceptual building 
blocks for the dissipative gadgets considered here are: an idea 
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FIG. 1. a) Depiction of an initialization gadget. The central circle 
depicts the target qubit, the outer circles depict the auxiliary qubits. 
The red lines connecting the timer qubit to the auxiliary qubits il- 
lustrate the conditional state preparation local Liouvillians £^ p , and 
the dotted circles around the auxiliary qubits illustrate unconditional 
amplitude damping Liouvillians Cf; i.e. £ ini = £ cp + £ ad . b) 
Depiction of the linear timer gadget. The circles represent the timer 
qubits, the lines connecting the timer qubits are the local Liouvillians 
£j Ut > 3 — 1> • • • 5 N — 1, and the boxes illustrate the initialization de- 
vices from a). 

for initialization, and a particular mixing behavior called the 
"cutoff phenomenon". Throughout, we will assume that the 
open system dynamics are properly modeled by a Markovian 
master equation with a local Liouvillian, so that 
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where each Lindblad operator Lj has a local support which is 
independent of the system size. 

The initialization gadget. For initialization to be success- 
ful, the qubit must (i) be prepared in the desired state in a time 
which scales at most polynomially in the gadget size M while 
the error decreases exponentially with M, and (ii) the state 
preparation must stop acting after a time which is specified by 
the desired dissipative process. We show how to initialize a 
single qubit in the |0) state, with a preparation time logarith- 
mic in the gadget size. For this, consider a star shaped graph 
as in Fig. 1 a), and label the central qubit as c and the auxiliary 
qubits as j = 1, . . . , M. The central qubit - the one we want 
to initialize - is surrounded by M auxiliary qubits. The Liou- 
villian of the gadget can be written £ im = £ ad + £ cp , where 
£ad re p resen ts a collection of amplitude damping channels 
acting independently on the M auxiliary qubits with rate lj, 

Lf = vM0fc)(l fc |, (1) 

for k — 1, . . . , M, and £ cp represents the preparation of the 
central qubit in state |0 C )(0 C |, at a rate T, conditioned on the 
auxiliary qubits being in the state 1 1&) (1& |: 

L c * = Vf|O c )(l c |0|l fc )(U|. (2) 



Intuitively, after a time tM = O (log M/lj), the M auxiliary 
qubits will all have relaxed to |0) to within exponentially small 
error, so that the state preparation essentially stops acting. We 
set out to show that for all but an exponentially small fraction 
of input states, the central qubit will be in the desired state |0 C ) 
after the time tM with an error that is exponentially small in 
M. Importantly, we point out that the dynamics of the initial- 
izer gadget acts exclusively on the diagonal elements of the 
system's density matrix. This implies that no off-diagonal el- 
ement contributes to the mixing time analysis, and hence we 
can assume that the initialization process is essentially classi- 
cal. 

Theorem 1 (Initialization) Let p be an arbitrary input state, 
and let C ini be the Liouvillian for the initialization gadget. 
If there exist <5, c > and a subset S C {1, . . . , M}, with 
\S\ = cM such that (lj\pj\lj) > S for all j G S. Then for 
any e > 0, there exists a r = O (log (M)) such that for all 

t > T, 

<lc|Pc(*)|lc) <Me^ M + £, 

where pj is the restriction of p to subsystem j, and p c (t) = 
tr aU xe^ (p) is the partial trace over the auxiliary qubits, 
and p is some positive constant which depends on {c, cj, T}. 

The precise form of p is given in the Appendix along with 
a proof. The theorem states that if a sufficient number of aux- 
iliary bits have a non-zero overlap with the |1) state, then the 
stationary state of the central qubit will be exponentially close 
to |0 C ). It also states that the system equilibrates in a time 
0(log (M)). It becomes clear in the proof that the assump- 
tions on the initial state can be further weakened, without af- 
fecting the result too severely. 

The construction of the initialization gadget can be ex- 
tended in several ways. The first convenient generalization 
is to initialize a given Lindblad operator rather than a state. 
This can be done simply by initializing a state, say |0) c and 
then conditionally attaching a Lindlbald operator L to |0) c . 
That way, L only starts acting when |0) c has been properly 
initialized. A second generalization is to extend the geometry 
of the initializer to graphs other than the star. In particular, if 
the initializer is constructed as a tree, with the intended ini- 
tialized state as the stem, then the initialization time can be 
engineered to be any polynomial in N. This is very useful 
when one wants to initialize a process which only mixes in a 
time polynomial in N. Furthermore, we point out that the gen- 
eral tree structure for the initialization gadget also allows us 
to embed it in a 3D lattice with nearest neighbor interactions 

ma. 

The timer gadget. We will now show how to trigger a dis- 
sipative operation at a particular point in time, and by con- 
catenation, how to perform sequential dissipative operations. 
Intuitively, one might be tempted to introduce timescales by 
simply imposing different decay rates 71 <C 72 <C • • • <C 7at 
before each desired dissipative operation Cj . But it becomes 
evident that either the convergence time or the error of such a 
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process will scale poorly in N. A way around this problem is 
to further exploit the "cutoff phenomenon" fT2l . A dissipative 
process is said to exhibit cutoff-type behavior |[T4ll . if there ex- 
ists an initial state which does not converge to stationarity for 
a long time, and converges exponentially fast (as a function of 
the number of qubits) after a specific time. 

We provide perhaps the simplest example of a timer gad- 
get: a line of N qubits such that the last one exhibits a cutoff 
from | ljv) to | Oat) in a time O (N). Each qubit of the timer 
must first be initialized in the state which will lead the chain 
to exhibit a cutoff. As depicted in Fig. 1 b), the timer gad- 
get consists of a main line of N "timer" qubits, connected 
by (N — 1) local Liouvillians C cut with Lindblad operators: 
Lf- = ^y|0 j ><0 J -|®|0 j+ i>(l J -+i|, forj = 1, . . . ,N-1. The 
local Liouvillian C^ ut acts as an amplitude damping channel 
with rate 7 on qubit j + 1 if qubit j is in 1 j ) and acts trivially 
otherwise. 

The initialization gadget guarantees that we start in the de- 
sired state vector |0 O ) •= |0) <8> |l)^ ) ( Ar - 1 ); i.e., the first qubit 
is in |0) and all of the other qubits are in |1). The qubit which 
will act as a trigger is the last one, so we want to estimate 
how long it will take before the N'th qubit is in |0), assum- 
ing proper initialization. This can be calculated explicitly. Let 
trjv_i be the partial trace over the first N — 1 qubits, then 
the following proposition proves cutoff behavior of the timer 
gadget in arguably the strongest possible sense. 

Theorem 2 (Timing) Let N be the number of timer qubits, 
then for t N (x) = (N + xy/N)/% 

i||tr 7 v_ 1 (e^^ cut ^o))-|0 A ,)(0 A ,||| 1 = l-*(s)+0(l/VAO, 

where §(x) = e~ f2 ^ 2 dt/ \/2ir is the cumulative normal 
distribution function. 

We note that for any real c > 0, it immediately follows 
from Theorem [2] that 

lim (0 N \trN-i[e ctN{0)CCU \fo)}\0 N ) = 0, fore < 1, 

TV— ^00 

lim (ON\trN-i[e ctN(0)CCU \ct> )}\0 N ) = 1, forc> 1, 

TV— ^OO 

which is the definition of a cutoff given in Ref. [14]. However, 
Theorem [2] makes a stronger statement yet. It gives us a cut- 
off time (tjsf = N/j), as well as the accurate cutoff window 
(At at = vN), and the deviation away from the cutoff time 
in the form of the cumulative distribution function <£. This 
last piece of information is especially useful as it allows us 
to estimate the inaccuracy of the timer, and analyze the error 
propagation for multiple uses of the gadget. 

It follows that if the state of the lattice is properly initial- 
ized to |0o), then an operation conditioned on the last qubit 
will only start taking place after a time of order t n = N/ 7, 
while it is exponentially suppressed before that. The trigger 
will operate on a time scale of order y/N, and the error in the 
triggering will be given by the tail of a normal distribution. In 
practice this means that a polynomial number of dissipative 



operations can be triggered each in linear time, with a linear 
time interval between them, while accurately triggering each 
operation with exponential precision. Rigorous error analysis 
is given in the Appendix. Finally, we note that the above lin- 
ear timer gadget displays three dissipative time scales, given 
by: 7 the rate decay of the timer Lindblad operators, F the 
decay rate of the initialization Lindblad operators, uo the de- 
cay rate of the amplitude damping operators on the auxiliary 
qubits of the initialization gadget. All of these time scales are 
constant in the size and in the number of the gadgets, so we 
can choose them in such a way that they clearly separate active 
dissipative processes. For optimal functionality we then want 
7 <C uo <C T, so that the timer Lindblad operators effectively 
do not start acting before 0o is initialized. We also emphasize 
that the timer gadget analyzed above is in no way unique, and 
several other models could be better suited to given situations 
lfT6l . We now turn to applications. 

Dissipative one-way quantum computation: We sketch how 
to perform D1WQC using the timer and initialization gadgets. 
We need to be able to perform the following three main steps: 

(i) initialize a two dimensional lattice in a cluster state [13], 

(ii) perform conditional projective measurements in a time or- 
dered fashion on the cluster state which simulate a universal 
gate set (here CNOT and arbitrary single qubit rotation) and 

(iii) propagate the classical boolean phase information from 
one measurement to another. All three of these steps have to 
be performed purely dissipatively, in such a way that the sta- 
tionary state is essentially unique, and that it is reached in a 
time which scales polynomially in the system size. It turns 
out that the obvious bottlenecks in the program - initialization 
of the cluster state and timing of operations - are solved by 
the timer and initialization gadgets introduced above. It can 
also be seen that, by construction, every Lindblad operator is 
geometrically local. 

Indeed, as shown in Ref. [14], a graph state lITTl of N qubits 
can be prepared dissipatively in a time which scales at most as 
O (log N)\ using 5-local Liouvillians in case of a 2D-lattice. 
Moreover, the state can be initialized by conditionally con- 
trolling each Lindblad operator of the graph state preparation 
with an initialization gadget. In particular, we can use the star 
geometry above, as the graph state can be prepared in time 
O (log N). 

Projective measurements can also be very naturally cast 
into a dissipative setting. Let {|£fc)} be the orthogonal basis 
associated with a given projective measurement on subsystem 
s, and let r denote the classical registry on which the mea- 
surement is stored, then the master equation describing the 
measurement is 

c meas {p) = £<akia)i s ® l&X&l - p- 

k 

Note in particular, that measurement is on the same footing as 
logical gates in this model, and hence, all resources are treated 
equally. 

In D1WQC it is also necessary to perform measurements 
conditional on the previous measurement outcomes. This can 
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be performed by keeping track of the pseudo-measurements 
on a classical information bus which is updated by condition- 
ing. This bus only needs to keep track of two bits per measure- 
ment and can be fed into future measurements by condition- 
ing. This can be performed locally on a three-rail bus. Here, 
each step in the D1WQC is triggered by a global clock (timer 
gadget). In Ref. fT3ll it was shown how to perform CNOT 
and one-qubit rotations in the 1WQC model. The CNOT 
gate can be performed by simultaneous measurements on nine 
qubits with two input and two output qubits. Arbitrary one 
qubit measurements can be performed on a line of five qubits, 
but the measurements have to be performed sequentially with 
feedback. We refer to Ref. fT3ll for a full description of these 
schemes. 

Rather than performing a full analysis of the dissipative ver- 
sions of these schemes, which is beyond the scope of this ar- 
ticle, we will give a detailed analysis of a subroutine, which 
already contains essentially all of the main ingredients of this 
model: qubit state transfer along a line. Consider a ID chain 
of odd-n spins, where we want to transfer the state of the 
first spin to the last spin by directed measurements. We as- 
sume that the initial state is already encoded in the cluster 
state at the input. In the standard setting this is done by first 
preparing the chain in the state \(pi n ) <S> |+) (g) ( n_1 \ where 
|d=) := (|0) ± \ l))/y/2, and then applying an entangling uni- 
tary operation Sj — |0j) (Oj | + |lj) (lj | <g) <rj +1 onto the state. 
The cluster state with \(fi n ) encoded in the input will be de- 
noted \cpi n )c '= S\cpi n ). We will not discuss how to initialize 
the input state dissipatively, but the generalization should be 
obvious. 

Now, in the standard state transfer by measurement, one 
would measure qubits 1 through n — 1 in the cr x basis and 
obtain measurement outcomes {si, . . . , s n -i} where Sj G 
{0,1}, and Sj = corresponds to the projection and 
Sj = 1 corresponds to the |— j) projection. The final state vec- 
tor after n—1 measurements is | out ) = |si) ® ■ • -® |s n -i) ® 
l^out)- It is not difficult to see that \(p ut) = U§\(fi n ), where 

TTSt S1+S3H \-S n — 2 S2~\ \-S n -l rpri 

= (J x cr z . The sum over measure- 

ment outcomes should be interpreted as being mod(2). In 
order to recover the original state one needs to perform the 
inverse unitary (U§)^ onto |yw)> or equivalently, if one is 
only interested in the measurement outcomes of \(f[ n ) after 
the state transfer, one can simply reinterpret the results based 
on the measurement outcomes {sj}. In order to illustrate the 
analogous dissipative construction, we show here how to per- 
form the transfer of a state along a three qubit chain. The 
full treatment of the n qubit state transfer is provided in the 
Appendix. 

Consider a chain of three qubits, labelled {1, 2, 3}, together 
with two auxiliary qubits, labelled {4, 5}. As in the stan- 
dard measurement based state transfer, we prepare the clus- 
ter state of three qubits with |(/? in ) encoded in the first; i.e. 
|<An)c- We will then act on the five qubits of the system with 
two maps sequentially; i.e. act with e tC until equilibrium is 

. n B 

reached, and then act with e 1 ^ . Both of these semigroups can 



be seen to mix rapidely, as they are conditional depolarizing 
semigroups with few conditioners. The sequential application 
of the maps can be performed by attaching each to a timer 
qubit which prescribes a specific starting time to each map. 
The Liouvillians C A and C B are given by the following Lind- 
blad operators, = V^I+jX+jl ® |0j+3>(lj+3|, L^j = 
(8) |lj + 3)(0j +3 |, where j = {1,2}. The second 
Liouvillian has Lindblad operators Lf = ^/uJa^ (g) |0, 0) (0, 1 1 , 
Li = V&Tf |0,0)(1,0|, Lf = V^rfaf |0,0)(1,1|, 
where the bras and kets with zeros and ones always refer to 
the auxiliary qubits. These maps will do the following: e tC 
takes \(pi n )c to an even superposition of 

| + , +) ® \<Pin) ® |0, 0), | + , -) <8> (<T z \<Pin)) ® |0, 1), 
|- +) ® Kl^in)) ® |1, 0), |- -) (<T x <T z \ipin)) ® |1, !)• 

In other words, it performs a measurement on qubits one and 
two but keeps a record of each outcome (on the classical reg- 
istry {4, 5}). The second map e tC acts as the unitary correc- 
tion U§ by rotating the superposition of states above onto the 
vector |0out)(0out| = li,2 ® |^in)(^in| |0,0)(0,0| , from 
which \ifin) can be read off directly. 

The CNOT and single qubit rotation are performed in much 
the same way with a slightly more involved classical bus for 
measurement and feedback (see the Appendix for the n-chain 
qubit state transfer). Likewise, concatenation of elementary 
gates can be handled efficiently, as the classical bus only ever 
has to keep two bits of information per updated qubit. 

Summary and outlook. In this work, we have introduced 
a basic toolbox for timing dissipative operations with a time- 
independent Liouvillian. We have shown that dissipative oper- 
ations can be performed sequentially with only a polynomial 
overhead in qubits or in Lindblad operators. We also show 
how to initialize (i.e., switch on and off) dissipative opera- 
tions, which then leads to full dissipative state initialization. 
We have provided the necessary error analysis to show that 
these gadgets can be applied sequentially while remaining ef- 
ficient. Finally, we outlined how these gadgets can be used to 
perform a dissipative version of one-way computation, hold- 
ing the promise of geometrically local dissipative computa- 
tion. The gadgets proposed here give rise to a versatile toolbox 
for constructing novel schemes of dissipative quantum infor- 
mation processing, opening up avenues for dissipative engi- 
neering. Finally, we mention a very promising applications 
for our gadgets: dissipative error correction, where the inter- 
play between active and passive operations become critical. 
A fully passive error correction scheme remains elusive, but 
clues suggest that dissipative engineering might be the right 
framework (see Appendix). It is the hope that the new per- 
spectives opened up here by introducing the possibility of se- 
quential schemes including measurement and processing trig- 
ger further research in this direction. 
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APPENDIX 

In this section, we prove Theorem 1, and we analyze an 
instance of a physically reasonable input state for the initial- 



ization gadget. To start with, we introduce some notation and 
state some basic observations which will be useful in the fol- 
lowing. First note that C lnl commutes with both the projection 
into computational basis and any permutation of the M outer 
qubits. So if we want to compute the overlap tr[Ae tC11Y \p)} 
w.r.t. any A which has the same symmetries, we can replace 
p by its symmetrization which has the form 

M 1 
m=0 o=0 

Here, the ' a are probabilities and p% = \a c ) (a c | ®Pk , where 
Pk is the normalized projection onto the set of all states of the 
outer M qubits with exactly k excitations, i.e., k qubits in |1). 
Without loss of generality, we can always restrict ourselves 
to the abelian algebra generated by the <^'s. In this way the 
problem becomes entirely classical. We will denote the con- 
ditional expectation, i.e., the projection onto this algebra, by 
V. 

Proof of Theorem 1 

Proof. We observe first that 

(l c |p c (t)|l c ) = 1 - (0 c |p c , D (t)|0 c ) < \\\p D (t) - ^||! 

where poif) = V o e tCini (p), p c ,D(t) is the corresponding 
reduced density matrix of the central qubit, and the stated in- 
equality uses monotonicity of the trace norm under the partial 
trace. Now, the stationary states of e tC ™ are of the form 

^®|0)(0|® M , 

where ^ = /3|O c )(O c | + (l-/3)|l c )(l c | and/3 e [0,1]. There- 
fore, consider the following two projections onto stationary 
states, 

T v (p) = (1 - UMtptiMMl + U\p\<Po)\<Po)(<Pol 
T 00 (p)= lim e tcin \p), 

t— >-oo 

and observe that o — T^. Then, using the triangle 
inequality for the trace norm, 

< \\(T v oe tciIli )(p D ) - ipX + \\(T v oe tciai - e t£ ' n ')(p D )||i. 

The second term can be bounded as 

K T v0e t^ _ e tc- ){pD)h 

< || (T v o e* £ini - T v o TooXpd) II ! + || (e* £ini - T^) (p D ) \\ 1 
<2||(e t£ini -T 0O )( /0l) )|| 1 <6Me-*- 

with the last inequality following from Lemma [4] 

The remaining term can be bounded by using Lemma[3j and 
writing pr> in terms of the basis • Given our assumption 
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that there exists a S > and a subset S C {1, . . . , M}, |5| = 
cM for some c > 0, such that (lj\pj\lj) > S for all j G 
5, we deduce that the probability of measuring all zeros in 
the input state is upper bounded by (1 — S) cM . Similarly, 
the probability of measuring k ones (k < cM) and all the 
rest zeros is upper bounded by (1 — S)( cM ~ k \ Thus, writing 

Pd = EfeioteM + VkVk)' and usin S Lemma [| with £ = 

cj/(7 + cj), we get 

M 

<^|e t£, %z>)l¥>o><£p& fe 

M(l-c) M 

< £ ^ fc + £ 

A:=0 fc=M(l-e) + l 

M(l-c) M 

< ^ (l-J) cM -^+ ^ ^ 
fc=0 fc=M(l-c) + l 

< Me^ M , 

where /i is the smallest (in absolute value) of all the exponents 
occurring in the sum. Clearly, every exponent is negative and 
proportional to M. Finally, we can write the whole expression 
as 

(UpS)\1c) <M(Ce-^ + e-^ M ). 

Choosing t > log (MC/e) /uj then completes the proof. □ 

One of the crucial properties of the initializer gadget, which 
is not transparent from Theorem 1, is that the gadget essen- 
tially stops acting on the central qubit after a time of order 
log M/oj. See the final comments for more on this issue. 

Lemma 3 Let ^q ' 1 ^ and C ini be defined as earlier, then 



<^|e* £ini (^)|^) = (l 



-t(uj+T)\k 



UJ 



Proof. To see this, consider the action of the initializer channel 
e tCini , where £ mi = £ ad +£ cp , on the pure product input state 
Lp\ . It is easy to see that 

£ ad (<^) = M^-i - <fi), £ a< V fe ) = MvLi - <pI). 

We want to know what the overlap with the stationary state 
ipl is after a time t. It is clear from the above equations that 
once the central qubit is in |0 C ), it never exits it. Since such 
states have zero overlap with p\ s , we can restrict ourselves to 
dynamics involving states with control bit |l c ). The reduced 
dynamics read 

C™{<p\) = -k{{T + w)ip\-w<p\_i). 

For simplicity, we assume that uj = T. The proof carries 
through almost unchanged when F ^ uj. In this case the dy- 
namics further reduce to 



£ ini (^) = -rfc(2^ 



Now, as we iterate £ im , we get an expansion in the basis {tp\}, 

m 

(>c ini ) m (^) = (-rr^^0> J 1 , 

where f^(j) are the coefficients of the expansion. It is clear 
by construction that f^(j) = if m + j < k, and that they 
obey the following recurrence relation 

/™+i(i) = -(j + + 1) + 2j/* (j). 

The coefficients f^(j) can be calculated explicitly, and are 
seen to be 



^(i) = ^E(z)(- 1 )'( 2 

1=0 v 



, +1 (?-l)!2i- 1 



, for j > 1, 



{I -3V- j! 

/-(°) = iE Cf){-i) l (ur, fori = 0. 



Now, we would like to know the overlap with (fl of the time 
evolved state ip\ . This is given by 

(^ie t£im (^)i^) = 2 ^|(£ !ni r(^)ko> 



m=0 



= E 



(-rt) m , fc 



ra=0 



m! 



I «/ m 



/£(o) 



m=0 



1=0 



1 

2^ 



(1 - e" 2 -) 



By proceeding in exactly the same manner, we get for uj ^ Y 
that 



; (l_ e -(r+a;)t ^fe 

cj + r 



□ 



Lemma 4 || (e^ - T^) {po)^ < 3Me~ 



Proof. We will regard any map as one acting on the abelian 
algebra in the range of V. If we define rj(T) := sup p \\T(p) — 
Too (p)\\i where the supremum is taken over all density matri- 
ces, then our task is to bound r/(e t£mi ). To this end note that 

we can write £ ini = Y,f=i &t where [4™' 4™] = if each 
Of 11 = + Cj d acts non-trivially only on the central and the 
j'th qubit. Since the Liouvillian is now a sum of commuting 
terms we can exploit Thm. 6 from Ref. ifTH which gives 

M 

i=i 

Cf l acts non-trivially only on two qubits, i.e., it is of the form 
Cf l = A (g) id, where A is a two-qubit Liouvillian. Now we 
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use that r] (e th (8) id) — r] (e £A ) . In order to compute the latter, 
we realize that e tA can be represented as a 4 x 4 stochastic 
matrix, for which r] can be computed using a computer algebra 
program which yields 



77(1 



'(2 



-tr 



u/(u + T)) < 3e~ 



□ 



Physically reasonable input state for the initialization gadget 

For completeness, we provide a wide class of physically 
reasonable input states which satisfy the conditions of Theo- 
rem 1, and explicitly give their scaling. Let p- in be the (classi- 
cal) input state, and write it as p- in = J2k=i Pk^k- We assume 
that pk follows a truncated discrete normal distribution with 
mean \i and variance a 2 . We further assume that there ex- 
ist strictly positive constants a, f3 > with a < 1 such that 
p = aM and a 2 = f3M. For the uniform input state 1/2 N , 
{pk} would indeed be a binomial distribution. We will now 
take the limit as N goes to infinity and work in the continuum. 
Hence we assume that {pk} obeys a truncated continuous nor- 
mal distribution with mean p and variance a 2 . The error in 
our calculation due to this assumption can be seen to decrease 
exponentially with N. We wish to evaluate 

lim (^|e t/:ini (p in )|^) = 

/l_ e -(r+-)^ fe 



k=0 



1 + T/lo 



Jo a 



where £ := 1 + T/cu and </>* r 6 is the truncated normal prob- 
ability distribution function restricted to the interval — 00 < 
a < b < 00 given by 



Here <p is the standard normal distribution function and is 
the cumulative normal distribution function. We show the fol- 
lowing: 

Lemma 5 Let a, f) > with a G (0, 1]. If p = aN and 
a 2 = f3N, then 



N 



t>o N ( ^ )Z x dx<Ce 



Nb 



for some N independent positive constants C, b. 

Proof We start by bounding the integral expression, and leave 
the denominator in the truncated normal distribution for later. 
Let z x := N(a - /31ogf) and z 2 := N(2a - /31og£). Then, 
a simple calculation shows that 



rN pN 



The integral on the right hand side is bounded above by 1, 
so if ^2 > 0, then the whole expression is bounded above 
by a function which is decreasing exponentially in N. There- 
fore, we can assume that z 2 < but this implies that z\ < 
because z\ < z 2 . Now, note that one can easily bound the 
complementary cumulative normal distribution from above by 
observing that for c > 0, c < x for any x G [c, 00], we get 



pOO pC 

/ e-* 2 / 2 dx< / 

J c J c 



~^l 2 dx 



-c 2 /2 



We use this to bound the integral as 



N 



dxq 



f x — [l 



/oo 



-e 2 £ . 



iV(a-/31ogO 

Hence, for any nonzero a, this expression decreases exponen- 
tially with N. The cumulative distribution functions in the de- 
nominator of the truncated normal distribution function can be 
bounded in the same way to show that the denominator is ex- 
ponentially close to 1. Thus putting all of the pieces together 
and bounding the pre-factor (which is strictly decreasing with 
TV), we get the desired result. □ 



Proof of TheoremU] 

Proof A simple calculation shows that 

N-2 , , k 

e t£Cut (^ ) = e"* + /(*)l°> • • • >0)(0, ■ • ■ , 0| 



k=0 



where fa is equal to the pure state with the first k + 1 qubits 
in the |0) state and the rest of the qubits in the |1) state, and 
f{t) is some function of t. Then, since tr[e tCCU (p)] = 1 for 
any p, and e tCCU is trace and separability preserving, we get 
that 



1, 



Jtrjv-ile"- ((l>o)]-\ON)(ON\\\i = l-(ON\tT N - 1 [e t ^ (fa)}\0 N 
and the expectation value can be evaluated as 

(0 N \tT N ^[e tcCU \fa)]\0 N ) = 1 - <l 7V |tr iV _ 1 [e^ CUt (0o)]|lAr) 

(fry)* 



N-2 ,, N i, 
k\ 

k=0 



1 



k\ 

T(N-l,t 7 ) 



T(N - 1) 
(N-2)\ ' 
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where 



r dr 



(3) 



is the Gamma function evaluated at y, T(y, s) 



r. 



e r dr corresponds to the upper incomplete Gamma 



function, and 7(2/, s) = / Q s r y 1 e r dr the lower incomplete 
Gamma function. The following relationships holds, 



7(j/,*)+r(y,5) =r(j/), 



(4) 



and T(n) = (n — 1)! when n is an integer. Now, it was shown 
inRefs. (231 El that 



7 (7V, TV + xVTV) = (N- l)\(&(x) + 0(1/ VN)). 

Thus, in the limit of large N, we get the desired result. □ 

Theorem[2]gives an estimate of the overlap with the incom- 
plete gamma function near x = as a function of N. In fact, 
the error term (0(1 /y/N)) is even exponentially suppressed 
in x. In particular, we treat the case when x is of order y/N 
below, where we consider successive triggering of dissipative 
operations. 



Imperfect initialization and error analysis 

The next natural question to ask is how much the timer is 
damaged by imperfect initializations? Given that the input 
state of the timer gadget is initialized with the initialization 
gadget, we can assume that it will be at least as close in trace 
norm to the desired input state 



0(s) = (|O)(O|(l-s)+|l)(l|s)®(|l)(l|(l-s)- 



for some worst case error e that was shown to scale expo- 
nentially in the number of qubits (e oc e _/iAr ) for some real 
fx > 0. We focus on a single worst case error e for simplicity 
of notation; it should be clear that a product input with differ- 
ent site-dependent error can be accommodated as well, which 
will not lead to significantly improved bounds, however. Let 
:= |0, 1, 1, . . . , 1)<0, 1, 1, ... , 1|, and let <j>^ denote <j> but 
with the bit in the k + 1 position flipped. Then, to first order 



m e, 



N-l 



k=0 

By the triangle inequality, and naively bounding the terms in- 
volved, 

(ONltiN^le^i^emON) « (0 N \tT N ^[e tccu \^}\0 N ) 

+Ne + (D(e 2 ). 



We now want to consider how the triggering error behaves 
under the concatenation of timer gadgets. For this we will 
invoke a variant of Theorem [2] What we want to show is 
that the concatenated triggering error is still negligibly small 
if we take N to be large. Assume perfect initialization, 
knowing from the previous section that imperfect initializa- 
tion only adds an error which is linear in e. Suppose we 
want to trigger L successive dissipative operations, where L is 
poly(iV). Each time step will be provided by TV linear timer 
gadgets. Then, for I = 1, . . . , L, we want to estimate the 
probability for operation I to be triggered in a time between 
h(l) = N(l - 1/2) /7 and t 2 (l) = N(l + 1/2) /j. As can be 
seen from the last step in the proof of Theorem 2, the proba- 
bility for correctly triggering an operation between t\(l) and 
^2(0 is one minus the probability of triggering before t\ (I) or 
after ^(0- This corresponds to 

7 (iVZ, N(l - 1/2)) F(Nl, N(l + 1/2)) 



r(Ni) 



r(m) 



These terms can be bounded as follows: It has been shown in 
Ref. EI that 

T(a,x) < 



x — a + 1 

Then, using Stirling's approximation, and dropping terms 
which are linear in N, we get 



T(Nl,N(l + 1/2)) 
T(Nl) 



< e- a ^ N po\y(N,l), 



where the degree of the polynomial is fixed, and a (I) is a 
strictly positive function of I. Thus, if we choose N polyno- 
mially larger than L, the probability of triggering an operation 
after ^(0 is exponentially small. Similarly, using the expan- 
sion 

/ \ a— 1 u 

TfrVPj = -x V_ 

T(a) ^ k\ ' 

v ' k=0 

we can bound the other term as 

7(NI,N{1- 1/2)) < e _ W)N 



T(Nl) 



poly(AM). 



For completeness, we give the functions a(l) and 0(1), 

a(l) = (1/2 - Z(log (I + 1/2) - log (l)))h, 
m = _(i/ 2 + j(log (I - 1/2) - log (Z)))/ 7 . 

Thus with N polynomially large in L, we can always guaran- 
tee that the probability for incorrectly triggering any of the L 
operations is exponentially small. 

Dissipative state transfer by measurement 



Hence, the behavior of the trigger is essentially unchanged as 
long as e is exponentially small in the system size. 



We describe here the general dissipative state transfer pro- 
tocol which will be slightly more involved than the three qubit 
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£? 



i i 

m 

3 = 1 



Bus qubits 



i I i 

m 



m 

i i i 



Measurement qubits 
Logical qubits 



and similarly for the second line of bus qubits. Finally, the 
recovery operation is almost identical to the three qubit chain: 

Lf = V^rJ-l <8> 10^-2, J -_ 2 >6<0 i _2, lj-2|6, 
Lf = y/utf_ i |0 j _2,0 j _2)6(lj-2,0 j _ 2 |6, 



FIG. 2. Depiction of dissipative state transfer of a qubit state along a 
chain of odd-n qubits. The lowest line of qubits represents the logical 
qubits prepared in a cluster state with the input state vector \(fi n ) 
encoded at the first site (see text). The red lines represent "dissipative 
measurements", whereas the red arrows depict dissipative classical 
information transfer along the bus, and (classical) dissipative binary 
logical operations from the auxiliary measurement qubits to the bus. 
The operations are performed from left to right sequentially in time 



chain case described in the main text, in that we need an extra 
register of qubit to keep track of the mod (2) sum of mea- 
surement outcomes corresponding to C/ff . The setup consists 
of odd-n logical qubits initialized in the cluster state vector 
|<An)c> represented in Fig. 2. Associated to the n logical 
qubits are n measurement qubits, and a classical information 
bus consisting of 2(n — 1) "bus qubits". The idea underlin- 
ing the n qubit chain state transfer is the same as for the three 
qubit chain in the main text, except that the recovery opera- 
tion will now be conditioned on the last two bus qubits only. 
This allows the recovery operation to remain local. Recall 
that, U§ — a2 X(J 2 z ' where j x — s\ + S2 + • • • + s n -2 and 
lz = $2 + ^4 + • • • + s n — 1 and addition is understood as 
being mod(2). The purpose of the classical information bus 
is to keep track of j x , j z as the state transfer moves forward. 

To write out the Lindblad operators, we label the logical 
qubits with a subscript I, the measurement qubits with a sub- 
script m and the bus qubits with a subscript &i ?2 . The Lind- 
blad operators performing the dissipative measurements are 
for j = {l,...,n- 1} 

L?j = M+M+Al®\0i)rn(lj\m, 



L 



3,J 



•j\m- 



The Lindblad operators which dissipatively update the bus 
qubits are for j = 2, . . . , n — 1 



for j odd, and 

for j even. Next, the information transfer from one state of the 
bus to the next is given by the following Lindblad operators 
for j = {2,...,n- 1}, 

^ = ^10^0^(0^,1^16!, 

L 2J = V^\0jAj)b 1 (lj:0 j \ bi: 



ua^a?^ <g> |0 j _ 2 ,0 j _ 2 )b(l j _ 2 ,l j _ 2 | b . 



The n — 1 measurements will project the system onto a su- 
perposition four sets of states. When one traces out the fi- 
nal logical qubit, the first element in the superposition will 
be \(pin), the second one cr z \(pi n ), the third one cr^l^in), and 
the fourth one cr x a z \(fi n }- The recovery operation then un- 
does the bit and phase flips to yield the desired final state. It 
is important to note that the bus update operations cannot all 
be performed simultaneously. So each step j must be per- 
formed successively. As the number of steps is polynomial, 
the scheme remains efficient. 



Dissipative quantum error correction 

In very much the same spirit as for D1WQC, we can con- 
sider dissipative error correction models by simulating the 
known active error correction models. A promising approach 
is to perform dissipative measurements on the qubits and pro- 
cess the information on a classical information bus in parallel 
before applying a corrective operation. The measurements, as 
well as the corrective operations, can be applied sequentially 
using a timer gadget, and the encoding and decoding can be 
performed efficiently with the initialization trick. In this way, 
it is tempting to think that one can implement Kitaev's toric 
code [20] dissipatively, hence yielding an efficient passive 
topological memory in 3D. However, a number of obstacles 
prevent us from certifying that this approach will surely work. 
First of all, despite a recent promising attempt | 21 ], there does 
not exist any protocol for actively correcting Kitaev's 2D toric 
code with fully local measurements and local classical infor- 
mation processing of polynomial depth. It is even question- 
able whether such an algorithm can exist in principle. A fur- 
ther obstacle standing in the way of a fault-tolerant dissipative 
quantum memory is that the gadgets which we introduced can 
also suffer errors. One way of overcoming this problem is by 
noting that the dynamics of the gadgets is essentially classi- 
cal, and hence can be protected using classical "dissipative" 
means. The prototypical example of such a classical dissipa- 
tive error correction is Toom's rule for asynchronous cellular 
automata, which was proved to scale favorably (22). Never- 
theless, this introduces a third drawback, which is that of in- 
terfacing the fault tolerant classical bits with the qubits; we do 
not know how to do this. Thus, although a number of obstacles 
stand in the way of a passive dissipative quantum memory us- 
ing the approach outlined above, we nevertheless believe that 
the flexibility provided by dissipative quantum processing, as 
illustrated through the timing gadgets, make it a very promis- 
ing framework for considering passive quantum memories. 
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Observations 

To conclude, we note two important facts related to the con- 
vergence of the schemes discussed here: 

(i) The timer and initialization gadgets which we have pre- 
sented and analyzed above have a constant gap, by con- 
struction. It turns out that in these processes, the second 
eigenvalue will be highly degenerate, with a degeneracy 
scaling linearly in N. This kind of spectrum is charac- 
teristic for processes exhibiting a cutoff phenomenon 

G2. 

This, and the application to D1WQC, implies that uni- 
versal quantum computing can be performed dissipa- 
tively with a global Liouvillian which has a constant 
gap in the system size. 

(ii) Once a state has been initialized, the initialization gad- 
get does not completely stop, but in fact keeps acting 
but with a strength which is exponentially suppressed in 
M. Hence, the action of the initializer after initializa- 
tion has taken place can be modeled by a quantum chan- 



nel given by f = ide + T(l -e), where e = G(e- rM ), 
and T is some noisy channel. Then, if we want to per- 
form operation St, we can bound the error by 

WfSt-StWoo^eWTW^Ke. 

Hence, the error will be exponentially small. An il- 
lustrative example of the above statement goes as fol- 
lows: Consider the rudimentary conditional operation 
acting on two qubits described by the Lindblad oper- 
ator L cond = v /7l°)( 1 l ® and assume that it is 
acting on the state p = (1 -e)\l, 0)(1, 0| 1|. 
In other words, the conditioning will be governed by e. 
It is easy to see that 

e tcc ° nd (p) = |1,0)(1,0| +e(|l,0)(l,0| + |0,1>(0,1|) 
+ ee-^(|l,l)(l,l|-|0,l)(0,l|). 

Thus, as long as e is very small, the conditional prepa- 
ration will be effectively stopped. In the case of the 
initialization gadget, this e is of order e~ tfl . 



